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ABSTRACT 

We employ X-ray stacking techniques to examine the contribution from X-ray undetected, mid-infrared- 
selected sources to the unresolved, hard (6-8 keV) cosmic X-ray background (CXB). We use the publicly 
available, 24//m Spitzer Space Telescope MIPS catalogs from the Great Observatories Origins Deep Survey 
(GOODS) - North and South fields, which are centered on the 2 Ms Chandra Deep Field-North and the 1 Ms 
Chandra Deep Field-South, to identify bright (S24fj.m > 80/iJy) mid-infrared sources that may be powered by 
heavily obscured AGNs. We measure a significant stacked X-ray signal in all of the X-ray bands examined, 
including, for the first time, a significant (3.2(t) 6-8 keV stacked X-ray signal from an X-ray undetected source 
population. We find that the X-ray-undetected MIPS sources make up about 2% (or less) of the total CXB below 
6 keV, but about 6% in the 6-8 keV band. The 0.5-8 keV stacked X-ray spectrum is consistent with a hard 
power-law (F = 1.44 ±0.07), with the spectrum hardening at higher X-ray energies. Our findings show that 
these bright MIPS sources do contain obscured AGNs, but are not the primary source of the unresolved 50% 
of 6-8 keV CXB. Our study rules out obscured, luminous QSOs as a significant source of the remaining 
unresolved CXB and suggests that it most likely arises from a large population of obscured, high-redshift 
(z> 1), Seyfert-luminosity AGNs. 

Subject headings: Galaxies: Active: Nuclei — Galaxies: Active: Infrared/X-ray — X-rays: Diffuse Back- 
ground 



1. INTRODUCTION 

Deep "blank-sky" surveys with the Chandra X-ray Obser- 
vatory (hereafter Chandra) have resolved ^ 50-90% of the 
hard (2-8 keV) cosmic X-ray background (CXB) into dis- 
crete sources (e.g., Bauer et al. 2004; De Luca & Molendi 
2004; Hickox & Markevitch 2006), with the resolved fraction 
decreasing with increasing X-ray energy (Worsley et al. 2004, 
2005). Subsequent spectroscopic observations have revealed 
that the majority of the sources are Active Galactic Nuclei 
(AGNs; e.g., Barger et al. 2003; Steffen et al. 2004; Szokoly 
et al. 2004). This is consistent with CXB synthesis models, 
which rely on the assumptions of the "unified" AGN model 
(see Antonucci 1993, for a review), that predict the observed 
power-law (F ~ 1 .4) shape of the CXB is created by the inte- 
grated X-ray emission from both soft (F ^ 1.8), unobscured 
and harder (F < 1.4), obscured^ AGNs (e.g., Comastri et al. 
1995; Gilh et al. 2007). 

Although deep X-ray surveys find the highest sky density 
of AGNs to date (- 7000 deg-^; Bauer et al. 2004), the de- 
creasing resolved fraction of the CXB with energy suggests 
that there exists an additional, highly-obscured AGN popula- 
tion that is missed in even the deepest X-ray surveys (Worsley 
et al. 2005). While the majority of the X-ray emission from 
these sources is attenuated, a small fraction of hard X-ray 
photons can penetrate the obscuring torus. In addition, hard 
(rest-frame > 10 keV) X-ray emission can be scattered into 
the observer's line-of-sight via Compton reflection, which 
comprises a relatively larger fraction of the total observed 
X-ray emission from a heavily obscured AGN. The fraction 
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of hard X-ray photons emitted from a heavily obscured source 
would likely be too small to identify the source individually, 
but the hard X-ray emission from many of these undetected, 
Compton-thick sources could comprise the unresolved CXB. 

Since the absorbed energy in obscured AGNs is re-emitted 
in the mid- and far-infrared, obscured AGNs should be bright 
Spitzer sources. The dust reprocessing of AGN accretion 
energy heats the dust to higher temperatures than can be 
achieved by heating via stellar processes, and this emission 
appears to be relatively isotropic (Lutz et al. 2004). The 
strong rest-frame 3-8/im thermal continuum produced by 
the AGN-heated dust can be used to separate luminous AGNs 
from starbursts using infrared color-color selection (Laurent 
et al. 2000; Lacy et al. 2004; Stern et al. 2005) and in- 
frared power-law selection (Alonso-Herrero et al. 2006; Don- 
ley et al. 2007). Once candidate AGNs have been identified. 
X-ray stacking analyses can be used to measure their average 
X-ray properties and determine their contribution to the CXB 
below the detection limit of individual X-ray sources. 

In this Letter, we use X-ray stacking techniques to exam- 
ine the contribution of MlPS-detected GOODS sources to the 
unresolved component of the CXB. In §2 we describe the 
X-ray and infrared data used in this study. Our analyses are 
presented in §3. Our discussion and conclusions are in §4. 
We use the cosmological parameters Hq = 70 km s~' Mpc~', 
f7M = 0.3, and f7A = 0.7. 

2. SAMPLE 

The 2 Ms Chandra Deep Field-North and 1 Ms Chandra 
Deep Field-South are currently the deepest 0.5 - 8 keV sur- 
veys (Giacconi et al. 2002; Alexander et al. 2003). The high 
spatial resolution of Chandra (^ 0."5 at the aimpoint) and 
Chandra's, low background make these images ideal for ex- 
amining the resolved fraction of the CXB. We use the X-ray 
catalogs of Alexander et al. (2003) for both the CDF-N and 
CDF-S, which contain 503 and 326 X-ray sources, respec- 
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Fig. 1. — Histogram of 24/im fluxes for the 2147 MlPS-selected sources 
in the GOODS-North & South fields (open histogram). The flux dis- 

tribution of MIPS sources with X-ray counterparts is shown for comparison 
(hatched histogram). 

lively, detected using wavdetect at a significance threshold 
of 10"^. We also consider the supplementary X-ray source 
catalogs for both fields, which add an additional 79 (42) 
sources to the CDF-N (CDF-S); these sources are detected 
with a looser threshold (10~^), but because they are coinci- 
dent with optically bright (R < 23) sources, the vast majority 
(> 90%) are likely real. 

To identify obscured AGN candidates, we use Spitzer data 
from the Great Observatories Origins Deep Survey (GOODS; 
Dickinson & GOODS Team 2005). The GOODS infrared cat- 
alogs consist of observations in the four IRAC bands (3.6, 
4.5, 5.8, and S.O^m; Fazio et al. 2004) and one MIPS band 
(24/im; Rieke et al. 2004). The contribution to the CXB from 
the IRAC-selected sources will be presented in a subsequent 
paper (Steffen et al. 2007, in prep). In this Letter we measure 
the contribution of the MIPS sources to the CXB. 

The GOODS MIPS observations are publicly available^ as 
pai-t of GOODS data releases DRl (GOODS-North) and DR3 
(GOODS-South). The GOODS-North and GOODS-South 
MlPS-selected catalogs contain 1199 and 948 sources, re- 
spectively, extending down to a 24/im flux limit of ~ 80/iJy. 
We cross-correlated the MIPS source positions with the afore- 
mentioned X-ray catalogs using a 3" matching radius. We 
found 253 (156) CDF-N (CDF-S) MIPS sources matched to 
an X-ray counterpart within 3" (with an expectation of ^ 20 
spurious matches). In Figure 1, we present histograms show- 
ing the distribution of 24/im fluxes for all MlPS-selected 
GOODS sources (open histogram) and for the MIPS sources 
with X-ray counterparts (hatched histogram). It is apparent 
that while the number of MIPS sources increases at fainter 
fluxes, the number of sources with X-ray counterparts in- 
creases slowly, so the fraction of MIPS sources with a de- 
tected X-ray counterpart decreases as one examines fainter 
24/im MIPS sources. This is consistent with the idea that lu- 
minous AGNs power the brightest 24/im sources, while the 
fainter sources can be powered by the more common dusty 
starbursts. 

3. ANALYSIS 

We used X-ray stacking techniques to examine the contri- 
bution of MlPS-selected AGN candidates to the unresolved 
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CXB. We calculated the total number of photons within a 3" 
diamater circular aperture, correcting for the fractional contri- 
bution from pixels only partially covered by the aperture. The 
local background is measured by extracting the total number 
of counts within a 30 pixel x 30 pixel (^ 15" x 15") box, cen- 
tered on a source and masking all of the known X-ray and 
MIPS source positions. 

To calculate the expected background in the desired aper- 
ture, the total number of background counts is scaled by the 
ratio of the summed exposure map values within the source 
aperture and the summed exposure map values in the local 
background region. This is similar to scaling by the relative 
sizes of the source aperture and background region, except 
this method accounts for variations in sensitivity due to chip 
gaps and edges, and HRMA mirror effects. Our simulations 
have shown that this method is both much faster and more 
accurate than measuring the local background using thou- 
sands of randomly placed apertures, as would be adopted in a 
Monte-Carlo approach. 

For each source, the encircled-energy fraction was calcu- 
lated using a tabular parameterization of the Chandra PSF, 
provided by the CXC, that gives the radius of a circular aper- 
ture for a given off-axis angle, photon energy, encircled- 
energy fraction, and azimuthal angle. ^ Using our desired 
aperture size, mean photon energy, and the sources' off-axis 
angles from the exposure-weighted mean aimpoints given by 
Alexander et al. (2003), we interpolated the HRMA FITS ta- 
ble to obtain the encircled-energy fraction for each source. 
Since the individual images used in the CDF-N and CDF-S 
mosaics were taken at a variety of roll angles, we removed 
the small azimuthal dependence by averaging the PSF param- 
eterization over all azimuthal angles. 

To calculate the total X-ray flux for all of the stacked 
sources, the total background-subtracted source counts (in 
photons) was divided by the sum of the mean exposure-map 
value within each extraction aperture (in cm^ s) to obtain 
the total photon flux for the stacked sources. This value 
was divided by the mean encircled-energy fraction, weighted 
by the mean exposure-map value for each source, and con- 
verted to energy flux by assuming a power-law spectrum 
with r = 1.4 and correcting for Galactic absorption using 
the X-ray opacity table of Morrison & McCammon (1983) 
and the exposure-weighted mean Galactic column density 
(Nhi = 1.16X 10^0 cm-2). 

In Table 1, we show the stacking results for the 638 
MIPS sources that lay within 6' of the exposure-weighted 
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Fig. 2. — (a) Stacked, adaptively smoothed, 4-6 keV X-ray thumbnail for 
the 638 MlPS-selected AGN candidates. Adaptive smoothing was performed 
using the csmooth algorithm in CIAO. The 3" diameter aperture used for 
source photometry is shown, (b) Stacked, adaptively smoothed 6-8 keV 
X-ray thumbnail. 



* Available at http://cxc.harvard.edu/cal/Hrma/psf/ECF/hrmaD1996-12- 
20hrci_ecf_N0002.fits. 



Constraints on the contribution of 24/im sources to the CXB 
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TABLE 1 

Results of X-ray Stacking Analysis of MIPS Sources 



Energy Band 


Total Counts 


Mean Background" 


Exposure Map 


Encircled Energy 


Area*^ 


Total Intensity 




%ofCXB'' 


S/N 


[keV] 


[counts] 


[counts] 


(mean) [cm^ s] 


Fraction'' 


[deg2] 


[ergs cm-2 s-' deg" 






l<y] 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 




[8] 


[9] 




Standard Chandra X-ray bands 


0.5-8.0 


6857.2 ±82.8 


4698.9 ±6.2 


3.86 X 10** 


0.556 


5.676 X 10-2 


(5.50 ±0.21) X 10- 


13 


2.35 ±0.09 


26.0 


0.5-2.0 


2904.6 ±53.9 


1429.8 ±3.4 


3.85 X 10* 


0.614 


5.930 X 10-2 


(1.64 ±0.06) X 10- 


13 


2. 16 ±0.08 


27.3 


2.0-8.0 


3952.6 ±62.9 


3269.1 ±5.2 


3.95 X 10* 


0.542 


5.559 X 10-2 


(4.22±0.34)x 10" 


13 


2.43 ±0.20 


10.8 


Narrow X-ray bands 


0.5-1.0 


1202.3 ±34.7 


679.2 ±2.4 


2.50 X 10* 


0.624 


5.950 X 10-2 


(0.60 ± 0.04) X 10- 


13 


1.99±0.13 


15.0 


1.0-2.0 


1702.4 ±41. 3 


750.5 ±2.5 


6.94 X 10* 


0.610 


5.920 X 10-2 


(0.85 ± 0.04) X 10- 


13 


1.86 ±0.09 


23.0 


2.0-4.0 


1491.7±38.6 


1088.5 ±3.0 


4.10 X 10* 


0.578 


5.826 X 10-2 


(1.36±0.13) X 10- 


13 


1.96±0.19 


10.4 


4.0-6.0 


1089.9 ±33.0 


927.8 ±2.8 


4.32 X 10* 


0.544 


5.588 X 10-2 


(0.97 ± 0.20) X 10- 


13 


1.73 ±0.36 


4.9 


6.0-8.0 


1371.0±37.0 


1252.3 ±3.2 


1.69 X 10* 


0.504 


5.463 X 10-2 


(2.83±0.89)x 10" 


13 


5.79 ±1.82 


3.2 



^The small error bars for the mean background result from the much larger area used to calculate this value, relative to the size of the extraction aperture. 
''This is the mean encircled-energy fraction weighted by the exposure map values of the sources. 

"This is the total area within 6' of the CDF-N and -S aimpoints, corrected for the masked areas within 2 X 90% encircled energy of the known X-ray sources. 
Assuming the CXB normalization of Hickox & Markevitch (2006). The CXB fraction is 11% higher assuming the normalization value measured by 
Revnivtsev et al. (2005) and 6% lower using the value from De Luca & Molendi (2004). 



mean aimpoints of the CDF-N or the CDF-S, and were out- 
side 2x the radius of the 90% encircled-energy fraction of 
the X-ray sources in the Alexander et al. (2003) main and 
supplementary catalogs to avoid contamination from known 
X-ray sources. The X-ray stacking was performed in the 
three standard Chandra bands, the full (0.5 - 8.0 keV), soft 
(0.5-2 keV), and hard (2-8 keV) bands, as well as five nar- 
rower, non-overlapping X-ray bands. 

From Table 1, the X-ray undetected MIPS source popula- 
tion makes up about 2% of the total CXB intensity below 
6 keV. The resolved fraction of the CXB increases to ~ 6% for 
the hardest X-ray energy band analyzed here. This is the first 
time that a statistically significant, stacked 6-8 keV signal 
has been found for an X-ray undetected population of sources. 
X-ray stacking analyses were used to measure the contribu- 
tion of the optical GOODS sources to the CXB, but no sig- 
nificant 6 — 8 keV signal was found (Worsley et al. 2006). To 
verify the authenticity of the 3.2(7 detection in the 6-8 keV 
energy band, we extracted thumbnail X-ray images centered 
on each X-ray undetected MIPS source and coadded them. 
In Figure 2, we show the resulting 15" x 15" stacked X-ray 
thumbnail for both the (a) 4 - 6 keV and (b) 6 - 8 keV bands 
with our 3" diameter aperture overlaid (black circle). It is 
clear from these smoothed images that there is a significant 
stacked X-ray signal within our extraction aperture for both 
X-ray bands {p = 6.9x 10~^ in the 6-8 keV band, assuming a 
one-tailed Gaussian distribution). Our 3" aperture was chosen 
to maximize the stacked signal within the aperture while min- 
imizing the background signal. While the peak of the stacked 
X-ray image is not precisely centered within our aperture, us- 
ing a larger aperture does not increase the measured S/N for 
the stacked 6-8 keV signal. 

While we detect a significant 6-8 keV X-ray signal for 
the stacked. X-ray-undetected MIPS population we do not 
know how this low-level X-ray flux is distributed among the 
sources. Is the stacked 6-8 keV signal dominated by a smaU 
number of bright 24^m sources, as one might suspect from 
Figure 1 ? To examine the importance of infrared flux on our 
X-ray stacking analysis, we break down our earlier stacking 
results by performing stacking as a function of limiting 24/im 
flux. In Figure 3, we show the resolved fraction of the CXB, 



the signal-to-noise (S/N), and the number of sources stacked 
as a function of limiting 24/im flux. Overall, it appears that 
the S/N of our stacked X-ray signal does gradually increase as 
we include the X-ray flux from fainter MIPS sources, which 
suggests that the stacked X-ray signal is not dominated by a 
small number of bright 24/im sources. 



4. DISCUSSION AND CONCLUSIONS 

We present in this Letter results from an X-ray stacking 
analysis of the X-ray-undetected 24/xm Spitzer MIPS sources 

MIPS Selection (6-6 keV) 




2.0 2.2 2.4 2.6 2.8 
L0G„(24 MICHON I'LUX [/iJy]) 



Fig. 3. — The total resolved 6-8 keV CXB fraction, the signal-to-noise of 
the stacked X-ray signal, and the total number of stacked sources as a function 
of 24 /im flux assuming the CXB normalization of Hickox & Markevitch 
(2006). Error bars (Icr) are shown for stacked sources that have S/N > 3, and 
3(7 upper limits are given for the stacked sources that are not significantly 
detected. 
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Fig. 4. — The redshift evolution of the observed-frame 2- 10 keV flux 
and 24fim flux density for three obscured AGNs, CX052 (Stern et al. 2002, 
dashed curve), NGC 1068 (solid curve), and the Circinus galaxy (dotted 
curve). The two dotted horizontal lines and one dotted vertical line denote 
the average flux limits for the 2 Ms CDF-N, 1 Ms CDF-S, and the GOODS 
MIPS catalogues, respectively. 

using the GOODS catalogs and the Chandra Deep Fields. 
We find, for the first time, a significant stacked X-ray sig- 
nal (3.2a) in the 6-8 keV band for an X-ray-undetected, 
AGN-candidate sample, suggesting that at least some of these 
sources harbor heavily obscured AGNs. Approximately 12% 
of the unresolved fraction (~ 50%; Worsley et al. 2005) of 
the 6-8 keV CXB can be attributed to these sources. The 
stacked 0.5-8 ke V spectrum has a power-law photon index of 
r = 1 .44 ±0.07, consistent with the slope of the CXB (Hickox 
& Markevitch 2006). The slope of the stacked X-ray spectrum 
hardens with increasing X-ray energy, suggesting that the flux 
from obscured AGNs dominates the stacked spectrum at high 
X-ray energies. 

We found evidence that there exist mid-lR-bright, heav- 
ily obscured AGNs that are not individually detected in the 
X-ray band, but we know Uttle about the general properties 
of these sources. By examining how the X-ray to 24/im 
flux relation for AGNs evolves with redshift, we can infer 
the AGN types that are detected in this mid-lR catalog. In 
Figure 4, we examine the evolution of the X-ray to 24yKm 
flux relation for obscured AGNs using the observed SEDs of 
three obscured AGNs: the luminous. Type 2 quasar CX052 
(distant Compton-thin, obscured QSO; L2-10 kev = 3.3 x lO'*^ 
ergs s~', intrinsic; Stern et al. 2002), NGC 1068 (extreme 
Compton-thick, luminous AGN; L2-10 kev ~ 10'*^ - lO"*"^ ergs 
s"\ intrinsic; Cappi et al. 2006); and the Circinus Galaxy 
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(Compton-thick, moderate-luminosity AGN; L2-i() keV ^ 10 
ergs s"', intrinsic; Matt et al. 1999). This is similar to Fig- 
ure 1 of Martinez-Sansigre et al. (2006), except here we uti- 
lize the observed SEDs of known obscured AGNs instead of 
using theoretical models. The rest-frame X-ray and IR fluxes 
were calculated by interpolating the SEDs provided by Stern 
et al. (2002) for CX052, and the NASA/IPAC Extragalactic 
Database (NED) for NGC 1068 and flie Circinus Galaxy. 

From Figure 4 it is apparent that a luminous, significantly 
obscured AGN such as CX052 would be detected in both the 
GOODS MIPS catalog and the deep X-ray surveys, even at 
z > 5, and would thus would not contribute to the unresolved 
portion of the CXB. Even if this source were Compton thick, 
it would easily be detected at 24/im and would be included 
in our X-ray stacking analysis. Given the smaU fraction of 
file CXB fliat we resolve at 6 - 8 keV, flie unresolved CXB 
at these energies cannot be attributed to a population of lu- 
minous, heavily obscured, "Type 2" Quasars at any redshift 
(unless Compton-thick QSOs are significantly fainter at 24/im 
than their less obscured counterparts). In addition, it is appar- 
ent that heavily obscured, low-luminosity Seyfert 2s without 
significant star-formation can be detected in the MIPS band 
at z < 0.8, but fall below flie GOODS detection threshold at 
higher redshifts. This suggests that the unresolved 6-8 keV 
CXB is not emanating primarily from a population of low- 
luminosity, low-redshift AGNs but could be from a popula- 
tion of low-luminosity, z > 0.8 AGNs. At these low 24/im 
fluxes it is difficult to separate AGNs from dusty starburst 
galaxies, which makes stacking analyses problematic due to 
the increased background signal in the 6 - 8 keV band from 
obscured starbursts. 

The average 0.5-2 keV flux for our MlPS-selected AGNs is 
^0.5-2 keV — 5.4 X 10"'^ ergs cm"^ s"', a factor of ~ 4.6 lower 
than the flux limit of the 2 Ms CDF-N (Alexander et al. 2003). 
Deeper X-ray observations will help to improve significantly 
the S/N of the stacked X-ray signal and will greatly improve 
our measurements of the 6-8 keV CXB from these sources. 
In addition, if it is possible to seperate better the obscured 
starbursts from the infrared-selected AGN candidates, the sig- 
nificance of the stacked X-ray detections should improve. We 
plan to address additional infrared AGN selection techniques 
in a subsequent paper (Steffen et al. 2007). 
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